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We revisit the question of the convergence of lattice perturbation theory for a pure SU (3) lattice 
gauge theory in 4 dimensions. Using a series for the average plaquette up to order 10 in the weak 
coupling parameter /3~\ we show that the analysis of the extrapolated ratio and the extrapolated 
slope suggests the possibility of a non-analytical power behavior of the form (1//3 — 1/5.7(1)) 10 ' 1 \ 
in agreement with another analysis based on the same asumption. This would imply that the third 
derivative of the free energy density diverges near (3 = 5.7. We show that the peak in the third 
derivative of the free energy present on 4 4 lattices disappears if the size of the lattice is increased 
isotropically up to a 10 4 lattice. On the other hand, on 4 x L 3 lattices, a jump in the third 
derivative persists when L increases, and follows closely the known values of f3 c for the first order 
finite temperature transition. We show that the apparent contradiction at zero temperature can be 
resolved by moving the singularity in the complex 1/B plane. If the imaginary part of the location of 
the singularity T is within the range 0.001 < T < 0.01, it is possible to limit the second derivative of P 
within an acceptable range without affecting drastically the behavior of the perturbative coefficients. 
We discuss the possibility of checking the existence of these complex singularities by using the strong 
coupling expansion or calculating the zeroes of the partition function. 

PACS numbers: 11.15.-q, 11. 15. Ha, 11.15.Me, 12.38.Cy 



I. INTRODUCTION 

It is widely accepted that the continuum limit of 
asymptotically free lattice gauge theories at zero temper- 
ature is obtained in the limit of arbitrarily small coupling. 
In practice, however, the relevant information about the 
continuum is obtained from a crossover region where both 
weak and strong coupling expansions break down, but 
where one can observe the onset of asymptotic scaling. In 
this region, observables such as the average of the elemen- 
tary plaquette are sensitive to the contributions of large 
field configurations. Such contributions can be modified 
by adding an adjoint term 0, 0, OJ 3 , a monopole chem- 
ical potential [a> El or by removing configurations with 
an action larger than some given value • These studies 
illustrate the fact that non-universal features or lattice 
artifacts seem generically present in the crossover region. 

On the other hand, one may hope that the universal 
features of the continuum limit may be obtained directly 
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from a weak coupling expansion. However, convergence 
issues need to be considered. The discontinuity in the 
plaquette average in the limits g 2 — + 0* precludes the 
existence of a regular perturbative series, and the decom- 
pactification of the gauge variables in lattice perturbation 
[i| should lead to asymptotic series 0- Despite these 
considerations, an analysis of the first 10 coefficients for 
the average plaquette P for the standard Wilson action 
in the fundamental representation 0, 0, 0, 0] suggests 
|14| a finite radius of convergence and a non-analytic be- 
havior of the form 

P^A^l/Pc-l/p) 1 -* . (1) 

In the present context, a ~ 0. If a > 0, the first deriva- 
tive of P diverges at (3 C and we say that the transition is 
second order since P is obtained by taking the derivative 
of the free energy with respect to (3. If — 1 < a < 0, the 
first derivative is bounded but the second derivative of 
P diverges and we call the transition third order. This 
type of singularity is not expected for the model consid- 
ered here (no adjoint term in the action). The divergence 
at C requires long range correlations and consequently 
a massless state. No such a state has been observed in 
glueball spectrum studies 0, 0, 0] . Consequently, a 
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natural strategy would be to try to falsify Eq. QJ. 

In this article, we test the validity of Eq. by di- 
rectly calculating the first and second derivative of P near 
the hypothetical C . We consider the minimal, unim- 
proved, lattice gauge model originally proposed by K. 
Wilson [18j. With standard notations, the lattice func- 
tional integral or partition function is 

Z = JJ f dUte-PXrV-WOMrrV,)) (2) 

with (3 — 2N/g 2 . Our study focuses on P and its deriva- 
tives which are defined more precisely in Sec. IIVI All 
numerical calculations are done with the fundamental 
representation of SU(3). 

In Sec. [n] we discuss various ways to estimate the 
unknown parameters in Eq. Q using the perturbative 
sereis. All the methods give results which are reason- 
ably consistent and show the robustness of the analysis 
of Ref. 0|- I n a ^ cases, a is very close to 0. If we use 
the scaling relation a = 2 — Dv with D = 4, a = goes 
with the mean field result v — 1/2. Note that a = is 
borderline between second and third order transition. In 
the mean field theory of spin models, the specific heat 
has a discontinuity. It is common to associate a discon- 
tinuity or a logarithmic divergence of the specific heat 
with a second order phase transition. At the other end, 
a discontinuity in the derivative of the specific heat is 
often called a third order phase transition. It has been 
observed in the context of the large- N limit of gauge the- 
ories in two dimensions 0], random surface models pfj| 
and spin glasses In 3D O(N) models, a ~ 0.11 for 
N = 1 (second order) and -0.12 for N — 3 (third order). 
In Sec. IIIII we discuss the possibilty of having asymp- 
totic series behaving temporarly as series with a finite 
radius of convergence and we discuss the effect of the 
tadpole improvement [22|. In Sec. IIVI we calculate the 
first two derivatives of P on L 4 lattices, and also on 4L 3 
lattices together with the average of the Polyakov loop. 
In Sec. [V] we show that the absence of peak increasing 
with the volume for the second derivative of P can be 
accomodated by moving the singularity in the complex 
1//3 plane. If the imaginary part of the location of the 
singularity is kept reasonably small but not too small, 
it is possible to limit the second derivative of P within 
an acceptable range without affecting drastically the be- 
havior of the perturbative coefficients. In Sec. IVII we 
discuss the possibility of checking the existence of these 
complex sin gula rities by using the strong coupling expan- 
sion j2j I24II25II or calculating the zeroes of the partition 
function |2al27ll2^|. 

II. SERIES ANALYSIS 

In this section we analyze the weak coupling series 
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FIG. 1: r m versus m for the three data set given in the text. 

We define r m = b m /b m _x, the ratio of two successive 
coefficients. The r m are displayed in Fig. ^ In this Fig. 
and the following, we compare three data sets. The first 
set is the 1995, order 8, series (empty circles) and 
the two other sets are the more recent [l^ order 10 series 
calculated on a 8 4 (stars) and 24 4 (filled circles) lattices. 
It is difficult to distinguish the ratios corresponding to 
the three data sets on Fig. ^ The first three coefficients 
are in good agreement with analytical results Jfjlllj . The 
apparent convergence of the ratios suggests [1J| that we 
consider a leading non-analytical behavior the form given 
in Eq. This will be our working hypothese for the 

rest of the section. 

One can estimate the unknown parameters |29j | with 
a fit r m = /3 C (1 + (a — 2)/m) For instance, the range 
3 to 10 for the ratios gives f3 c ~ 5.66 and a ~ 0.24 
while the range 4 to 10 gives (3 C ~ 5.69 and a ~ 0.18. 
The 1 /to 2 corrections were reduced with a modified form 
r m = f3 c (l + (a - 2)/(m + s)) in Ref. with a value 
s = 0.44. Using the range 4 to 10 again, one obtains 
r . — 5.78 and a ~ —0.01. The estimates quoted in Ref. 
[3 are (3 C ~ 6 x 0.961 ~ 5.77 and a ~ 0.01. 

Another way to reduce the uncertainties of order 1 /to 
in the estimation of /3 C is to use [3(j the extrapolated 
ratio (R m ) defined as 

R m = mr m - (to - l)r TO _! . (4) 

The values of R m are shown in Fig. |3 As the calculation 
involves differences of successive ratios mutiplied by the 
order, the result is sensitive to statistical errors. Fig. |3 
illustrates the fact that the 2000 series has much smaller 
statistical errors than the 1995 series. Fig. [21 also indi- 
cates that as the volume increases from 8 4 to 24 4 , R m 
increases by approximately 0.03 which is slightly larger 
than the changes in R m with to < 8. Our best estimate 
of j3 c is Rio — 5.74 for the 24 4 lattice. Given the other es- 
timates and the size of the volume effects, it is reasonable 
to conclude that (3 C = 5.7(1). 

The corrections of order (1/to) 2 in r m introduce an 
uncertainty of order 1/to in the estimation of a. It is 
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FIG. 2: Rm versus m for the three data sets described in the 
text. 
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FIG. 3: The normalized slope S„. 

possible to eliminate this effect by using [3(j the extrap- 
olated slope (S m ). For this purpose, we first introduce 
the normalized slope S m defined a 

S m = -m(m-l)(r m -r m _i)/(mr m -(m-l)r m _i) (5) 

The values of S m are shown in Fig. [3] again showing a 
much better stability for the 2000 series. The volume 
effects are smaller than for R m . They are typically of 
order 0.01 or smaller. It is known [3(| that S m ~ a — 
2 + K/m. Using a fit of this form, we obtain a = —0.03 
for the m = 4 to 10 values. One can then remove the 
1/m corrections by using |3(j the extrapolated slope S m 
defined as 

S m = rnS m - (m - l)S m -i ■ (6) 

It provides an estimator with corrections of the form 

S m = a - 2 - Bm- A + 0(m- 2 ) , (7) 

where A is an exponent corresponding to possible nona- 
lytical corrections to the leading term Eq. ^ In 3 dimen- 
sional scalar models A is related to irrelevant directions. 



FIG. 4: S m versus m. 

The values of S m are shown in Fig. 0J Again, the fluc- 
tuations are larger than for the unextrapolated quantity 
and the 2000 data is much more stable. For the 2000 
data on a 24 4 lattice, we have S9 = —1.91 S10 ~ —2.08. 
Based on other estimates, we conclude that |a| < 0.1. 

In summary, the analysis of ratios suggests that P has 
a non- analytic, power- like singularity 

P n .- a . oc (1/5.7(1) - l//?) 1 - ^ . (8) 

III. REMARKS 

It is not difficult to find an asymptotic series (which 
has a zero radius of convergence) for which the begin- 
ning coefficients suggest a power singularity. A simple 
example is 

poo 

Q{p)= / dte-^l-i/VCB/?)] 1 -" (9) 
Jo 

the ratios of coefficients of can be calculated exactly: 

r m = ((B + m)/B)p c {\ + {a - 2)/m). 

For B sufficiently large and m « B, we have 

r m ~ /3 C (1 + (a - 2)/m), 

which corresponds approximately to the situation ob- 
served for the series discussed above. On the other hand, 
for m >> B we have r m cx m whichs shows that for 
sufficiently large m, the coefficients grow factorially. 

It is also well-known that the convergence of pertur- 
bative series can be improved by using the tadpole im- 
provement [22J. One defines a new series 

K K 
rn— m— 

with a new expansion parameter 

1 Z^m=0 Um f- J 
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FIG. 5: First derivative of P versus (3. 



The new coefficients are shown in Table I. If we exclude 
the last coefficient which may not be very reliable, we find 
that series is alternate and has much smaller coefficients. 
The first six ratios are of order -1, which suggests a sin- 
gularity at negative (3r. A discontinuity in P at (3 ~ —22 
was observed in 8], however we have no interpretation 
for p R < 0. 

IV. DIRECT SEARCH FOR A SINGULARITY 

If we take Eq. JHJ at face value, it implies that the 
third derivative of the free energy density should have a 
singularity with an exponent close to -1. More precisely, 
we would expect: 

d 2 P/df3 2 oc (1/5.7(1) - l/P)- 1 - ^ . (12) 

For symmetric finite lattices with L D sites and periodic 
boundary conditions, the number of plaquettes is 

J\f p = L D D(D - l)/2 . (13) 

Using the free energy density 

/ = -(l/AginZ, (14) 

we define the average plaquette 

P = df/d/3 = (1/Ag (£) (15) 

with 

£ = 2(1 - (l/N)ReTr(U p )) . (16) 
v 

We also define the higher moments: 

-dP/dp={l/N p )[(Y?)-^) 2 ] , (17) 

and 

d 2 P/dp 2 = (1/Ag[(£ 3 ) - 3 <£) <E 2 ) + (S) 3 ] (18) 



It should be noted that there is some loss of precision 
in the calculation of the higher moments. For instance, 
in —dPd[3, the two terms are of order Af p but their dif- 
ference is of order 1. For d 2 P/d/3 2 , the three terms are 
of order Af 2 while their combination is of order 1. For a 
symmetric lattice with 10 4 sites, for instance, the second 
derivative of the plaquette will appear in the ninth signif- 
icant digit and the use of double precision is crucial. To 
give an idea, 400,000 reasonably uncorrelated configura- 
tions were necessary in order to get a signal significantly 
larger than the statistical fluctuations for the four points 
on a 10 4 lattice in Fig. For (3 > 5.7, the second deriva- 
tive seems to be of the same order or smaller than the 
fluctuations for the same volume. 

We have calculated the plaquette and its first two 
derivatives on L A lattices, with L— 4, 6, 8 and 10. The 
first derivative is shown in Fig. When going from 
L = 4 to L — 6, the peak moves right and its height 
diminishes. This situation is very close to Fig. 1 of Ref. 
[3l| in the case of SU(2). When increasing L to 8, the 
values around 5.7 slighly drops and then stabilizes for 
L = 10. 

The second derivative of the plaquette is shown in Fig. 
For comparison, we have also reproduced the first 
derivative over the same range and with the same scale. 
In general, we found a reasonable agreement between the 
values obtained by subtracted averages or by taking the 
a numerical derivative as found in |31| . The figure makes 
clear that peak seems to disapear to a level close to our 
statistical fluctuations when the volume increases. 

We have also calculated the same quantities on 4L 3 lat- 
tices with L = 4,6 and 8. The results are shown in Figs. 
00 and El The sudden jump in the second derivative of 
P persists when L increases. The location of the jump 
coincides with the onset of the average of the Polyakov 
loop. The critical values of (3 corresponding to a first or- 
der finite temperature transition are well-known for 4L 3 
lattices [53, H3, HH- For L sufficiently large, it occurs 
very close to (3 — 5.69. The locations of the jumps seen in 
the second derivative in Figs. 7-9 follow very closely the 
quantity the real part of the leading zero of the par- 
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FIG. 6: First and second derivative of P versus (3. 



tition function, of Table 9 of Ref. [23 (5.552 for L = 4, 
5.650 for L = 6 and 5.674 for L = 8). 

The fact that for a 4i 3 lattice, (3 C ~ 5.69 is very close 
to /3 C ~ 5.74 obtained with the extrapolated ratio from 
the perturbative series seems to be a coincidence. Had we 
taken N t to be 12, we would have obtained |3i| (3 C — 6.33. 
Note that on L 4 lattices, the onset of the Polyakov loop 
occurs at larger values of (3. For instance, on a 8 4 lattice, 
it is 0.1 near (3 = 6.6. However, we found nothing but a 
smooth decrease for the first derivative of P in this range. 



V. RESOLUTION OF THE PARADOX 

It is possible to regularize the singularity that, accord- 
ing to Eq. JU, would appear in d 2 P/d(3 2 by replacing 
the singularity at 1//3 C by a pair of complex conjugated 
singularities in the 1/(3 complex plane. If we denote the 
imaginary part of the location of the singularity by T, 
then some apparently difficult compromise needs to be 
achieved. On one hand, if Y is too large, it produces pe- 
riodic modulations of the ratios. Fig. Ogives no hint of 
such an effect. On the other hand, if T is too small it 
does not fulfill its purpose. The question is thus if wc 
can find a window of acceptable values for V. 

A simple alternative to Eq. (JIJ can be designed by as- 
suming that the critical point in the fundamental-adjoint 
plane has mean field exponents 2] and in particular 
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FIG. 7: The Polyakov loop, P, first and second derivative of 
P versus f3 for a 4 4 lattice. 



a = 0. We will further assume an approximate loga- 
rithmic behavior 



-dP/d/3 oc ln((l//3 m - 1//3) 2 



(19) 



on the axis where the adjoint term of the action is zero 
(the range of parameters considered here). 1/ (3 m denotes 
the value where the argument of the logarithm is maximal 
on this axis. This implies the approximate form 



d 2 P/d(3 2 



-C 



(l/(3 m - 1/(3) 



(3*((l/(3 m - l/(3f 



T 2 ) 



(20) 



The (3 3 at the denominator ensures that the series starts 
at (3~ 3 . The three unknown parameters can be approx- 
imately fixed by minimizing some linear combination 
(with positive "weights" ) of the square of the differences 
between the coefficients of the series Eq. (I20f) and the 
actual ones. We expect the coefficients of larger order to 
carry more information about the non-analytic behavior, 
consequently, more weight should be given to the large 
order coefficients. For instance if we minimize a \ 2 which 
is the sum of the square of the relative errors for the last 



6 



H 0.1 



0.0 
0.51 



0.45 



0.42 
0.30 



$ 0.20 

0.15 
2 



X 



0.48 



0.45 



0.30 ■ 

0.25 

0.20 
3 
2 



5 5.6 5.7 5. 

P 



— •— 4x6 3 




*V — •— 4x8 3 
\ 

* 

\ 


• 4x8 3 


\ 

\ 

• 

\ 
\ 


— •— 4x8 3 





5 5.6 5.7 5. 



FIG. 8: The Polyakov loop, P, first and second derivative of 
P versus (3 for a 4 x 6 13 lattice. 



FIG. 9: The Polyakov loop, P, first and second derivative of 
P versus f3 for a 4 x 8 3 lattice . 



four coefficients, we obtain (3 m ~ 5.78, T ~ 0.0058, and 
C ~ 0.15. Integrating twice, we obtain the approximate 
series with coefficients which can be compared with the 
c m given in Table II: 



P 



~ 0.44/T 1 
f32.27/T 5 
+2533/T 8 



0.85/T 2 

132.3/3" e 

11590/3" 



2.44/3" 3 + 8. 
f 568.9/T 7 
+ 54160/3" 10 



42/T 



(21) 



The agreement increases with the order. The relative 
error on the fifth coefficient is about 6 percent and keep 
decreasing to a level smaller than the numerical errors 
as the order increases. We have tried with many other 
weights involving the last 6 coefficients and found very 
little variations in the estimation of C (typically \6C\ ~ 
0.01) and /3 m (typically \S/3 m \ - 0.02). On the other 
hand, T varies more rapidly under changes of the weights 
in the \ 2 function. We found values of T between 0.003 
and 0.007. 

The stability of C and (3 m can be used to set a 
lower bound on T. Given that the approximate form 
of d 2 P/d[3 2 in Eq. ignj has extrema at 1/(3 = l//3 m ±T. 
As we do not observe values larger than 0.3 near (3 — 5.75 



(see Fig. |SJ| we get the approximate bound 
C 



< 0.3 



(22) 



This implies the lower bound T > 0.001. On the other 
hand, large values of T are also excluded. As we never 
found estimate close to 0.01, we conclude that 



0.001 < T < 0.01 



(23) 



We also performed calculations with an assumption 
similar to Eq. but with ((l//3 m - l/{3) 2 + r 2 ) 1+ ( Q / 2 ) 
at the denominator, for small positive and negative val- 
ues of a. We found very similar ranges of values for the 
unknown parameters and we were able to draw very sim- 
ilar conclusions as for a = 0. 

A puzzling aspect of Fig. [S] is that the maximum of 
the first and second derivatives are not located near 5.78 
but near lower values (5.55 and 5.63 repectively) . Using 
the parameters obtained with the procedure described 
above, we found that a typical value for the maximum of 
the second derivative of P is 0.1 or below. This is clearly 
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below our numerical resolution and below the maximum 
value 0.6 in Fig. This can be explained from the fact 
that the perturbative series becomes a poor estimate of P 
when j3 becomes too small. Using the parametrization of 
this difference (the " non-perturbative part of the plaque- 
tte") of Ref. and taking two derivatives, we estimate 
that the nonperturbative contribution to d 2 P/dp 2 can 

be approximately written as 2.84 x 10 11 x exp(— 16 J 3 13 ) 
for 5.6 < (3 < 5.8. This function takes the values 0.25 
at — 5.8 and dominates the perturbative part. It also 
takes the value 0.65 at p = 5.6 which is consistent with 
our numerical calculation. As the parametric form of the 
nonperturbative contribution is still being debated, this 
is not the last word on the question, however, it seems 
clear that the nonperturbative part plays a major role in 
explaining Fig. 



VI. CONSISTENCY CHECKS 

It would be interesting to check if the complex singu- 
larities suggested by the present analysis could be seen 
using independent methods. In this section, we discuss 
the possibility of using the strong coupling expansion and 
the zeroes of the partition function for this purpose. The 
strong coupling expansion of the free energy has been cal- 
culated up to order 16 in Ref. [H for SU{2) and SU(3). 
With a suitable rescaling of P, these series can be used 
to calculate P and its derivatives in power of P (with 
the normalization of Eq. @). In the case of SU(2), the 
complex singularities of the specific heat were discussed 
in various sets of complex variables that have a simple 
interpretation for small values of P [24l |25| and com- 
pared with zeroes of the partition function on a 4 4 lattice 
[26j. These zeroes were calculated using a re- weighting 
of Monte Carlo data at real /?. 

In order to provide a comparison, we will start with 
the case of SU(2) which is better understood. We con- 
structed the expansion in powers of P of P using Ref. 
|23[ . This serie is 1 plus an odd series in p. This can be 
seen as a consequence of the identity P(P) + P{~P) = 2 
derived in Ref. 0. The sign of the nonzero coefficients 
alternates. The ratio of successive odd coefficients seems 
to converge toward a value close to —0.4. This indicates 
a singularity at P 2 ~ —2.5, or in other words, a pair of 
purely imaginary singularities at P ~ ±il.6. We found 
36 Pade approximants of P with their singularity closest 
to the origin located nearby (within 0.1 in both coordi- 
nates) this estimate. More sophisticated estimators such 
as the ones of Sec. |n]do not improve the accuracy of the 
estimation. In [24J, the single loop variable J was used 
instead of p. Using a ratio analysis of the series in J, 
these authors concluded that there is a singularity near 
J 2 ~ —1/6. This corresponds to P ~ ±il.5, in good 
agreement with the estimate we gave above. 

We have inspected the poles of the Pade approximants 
for P and its first two derivatives in a horizontal strip 
defined by the conditions ReP > 1.5 and |7m/3| < 1. We 



found single poles on the real axis and pairs of conjugated 
poles rather far from the real axis. To fix the ideas, for 
the second derivative of P, the pair closest to the real axis 
in this strip is 1.79 ± i0.43 for a [3/9] approximant. Also, 
the poles on the real axis tend to cluster. For instance, for 
the first derivative of P there are 24 approximants with a 
pole between 2.1 and 2.5, and for the second derivative, 
13 approximants with a pole between 1.7 and 1.8. 

In Ref. [24| , an attempt was made to find evidence for 
complex singularities in a new variable z = 6 J 2 / (6 J 2 + l). 
It was concluded that the existing series was too short to 
see a clear departure from a singularity on the real axis 
(near P = 2.2). On the other hand, complex zeroes of 
the partition function were found |2a | on 4 4 lattice near 
P = 2.2 + ±i0.15. It is plausible that the imaginary part 
increases with the volume and possibly reaches the es- 
timate P ~ 2.2 + ±i0.3 based on a more sophisticated 
treatment of the high-temperature expansion |25|. We 
are not aware of any numerical check of this statement 
on larger lattices. However, the fact that the two inde- 
pendent estimates are close and that one expects the spe- 
cific heat peak to broaden with the volume, makes the 
existence of a singularity at P — 2.2(1) ±«0. 2(1) quite 
plausible. Consequently, it seems that the Pade approx- 
imants often provide a good estimation for the location 
of the singularities close to the origin but not for the ex- 
pected singularities close to the real axis that are farther 
from the origin. 

For SU(3), we analyzed the /? expansion of P following 
a similar procedure. With the exception of the third 
coefficient, the coefficients are nonzero. The signs do 
not show any obvious periodic pattern. If we plot the 
logarithm of the absolute value of the coefficients, we see 
that they approximately fall on a line, however there are 
exceptions. We already mentioned the vanishing third 
coefficient, in addition, the 11th (15th) coefficient falls 
significantly below (above) the linear fit. The slope of 
the linear fit is -1.49 if the coefficients 4 to 15 are used. 
This indicates a radius of convergence of approximately 
4.5. Due to the significant discrepancies from the linear 
behavior, the ratio analysis does not reveal any obvious 
information. 

The poles of the Pade approximants are quite irregu- 
lar. However, many approximants have a pair of complex 
conjugated poles near 5e ±ITr / 3 . We have found 12 Pade 
approximants of P for which the poles closest to the ori- 
gin were located close to these two points (within 0.3 
in each coordinate). An inspection of the poles of the 
Pade approximants for P and its first two derivatives in 
a horizontal strip defined by the conditions ReP > 4 and 
|7ro/3| < 2 shows single poles on the real axis and pairs of 
conjugated poles far from the real axis. For instance, for 
the second derivative of P, the pair closest to the real axis 
in this strip is 4.74±il.72 for a [5/5] approximant. Con- 
sidering that in the case of SU (2) we were unable to find 
the expected singularity near P ~ 2.2 + ±i0.2, it seems 
plausible that a similar phenomenon occurs for SU{2>). It 
would be interesting to repeat the above analysis using 
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changes of variables as in the SU(2) case. 

In Sec. IIVI we have already mentioned the fact that 
the leading zeroes of the partition function for SU (3) on a 
4L 3 lattice had been calculated |53 • On a 4 4 lattice, they 
are located near j3 — 5.55±i0.12. On larger symmetric L 4 
lattices, we expect that the imaginary part will increase 
with L. We are not aware of numerical calculations for 
L > 4. With the parametrization of Sec. |V] we predict 
complex singularities at j3 = f3 m ± «/3„r. For (3 m = 
5.75 and 0.001 < T < 0.01, we predict an imaginary 
part between 0.03 and 0.3. It seems feasible to check 
this prediction for 8 4 or larger lattices, using the method 
documented in Ref. |27| . 

It should also be noted that new methods have been 
developed to determine the order of a phase transition 
using the distribution of zeroes of the partition function 
and their impact angle in the reduced variables complex 
plane [2fij. It would be quite interesting to apply these 
methods to interpret the variations of zero distribution 
when an adjoint term is added to the action. 

VII. CONCLUSIONS 

We have shown that the apparent conflict between an 
hypothetical singularity in the second derivative of P sug- 
gested by the perturbative series and the absence of evi- 
dence for a peak with height increasing with the volume 
on isotropic lattices, can be resolved by moving the sin- 
gularity in the complex 1/ ' (3 plane. If the imaginary part 
of the location of the singularity T is within the range 
0.001 < r < 0.01, it is possible to limit the second deriva- 
tive of P within an acceptable range without affecting 



drastically the behavior of the perturbative coefficients. 
This picture seems consistent with a small value of a but 
our numerical analysis does not single out the mean field 
value a = 0. 

It would be interesting to calculate the zeroes of the 
partition function for SU(3) on 8 4 or larg er symmetric 
lattices using the the re-weighting [2(| l27j method. We 
expect that as the volume increases, the locations of the 
leading zeroes should stabilize at values 5.7(1) ± i0.2(l). 
Changes of variables analog to those used in 51/(2) j2J, 
|25| should be used in order to see if there is agreement 
with the other methods. It would also be interesting to 
follow the distribution of zeroes and their impact angle 
when an adjoint term is added using the results of Ref. 
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